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Abstract 

We consider bosonic atoms with a repulsive contact interaction in a trap potential for a Bose- 
Einstein condensation (BEC) and additionally include a random potential. The ensemble averages for 
two models of static (I) and dynamic (II) disorder are performed and investigated in parallel. The 
bosonic many body systems of the two disorder models are represented by coherent state path integrals 
Zi[J'], Zii\J\ on the Keldysh time contour which allow exact ensemble averages for zero and finite 
temperatures. These ensemble averages of coherent state path integrals therefore present alternatives to 
replica field theories or super-symmetric averaging techniques. Hubbard-Stratonovich transformations 
(HST) lead to two corresponding self-energies for the hermitian repulsive interaction and for the non- 
hermitian disorder-interaction . The self-energy of the repulsive interaction is absorbed by a shift into 
the disorder-self-energy which comprises as an element of a larger symplectic Lie algebra sp(4M) the 
self-energy of the repulsive interaction as a subalgebra (which is equivalent to the direct product of M 
times sp(2); 'M' is the number of discrete time intervals of the disorder-self-energy in the generating 
function Zi\J\). After removal of the remaining Gaussian integral for the self-energy of the repulsive 
interaction, the first order variations of the coherent state path integrals Zi\J\, Zii\J\ result in the 
exact mean field or saddle point equations, solely depending on the disorder-self-energy matrix. These 
equations can be solved by continued fractions and are reminiscent to the 'Nambu-Gorkov' Green 
function formalism in superconductivity because anomalous terms or pair condensates of the bosonic 
atoms are also included into the selfenergies. The derived mean field equations of the models with 
static (I) and dynamic (II) disorder are particularly applicable for BEC in d = 3 spatial dimensions 
because of the singularity of the density of states at vanishing wavevector. However, one usually starts 
out from restricted applicability of the mean field approach for d = 2\ therefore, it is also pointed out 
that one should consider different HST's in d = 2 spatial dimensions with the block diagonal densities 
as 'hinge' functions and that one has to introduce a coset decomposition 5p(4M)\?7(2M) into densities 
and anomalous terms of the total disorder-self-energy sp(4M) for deriving a nonlinear sigma model. 
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1 Introduction 

The original Bose-Einstein condensation of atoms in traps has been extended and performed in various 
manners yy-|6j. In this article we investigate a BE-system of bosonic atoms with repulsive interactions in 
a trap potential and include a random potential which represents a model of ensemble averaged disorder. 
We concentrate on the averaging procedure, derive the exact saddle point equation and describe continued 
fractions as solutions of the mean field equations which also comprise anomalous terms or pair condensates 
in the coset parts of the selfenergies. The exact saddle point equation is obtained from first order variation 
of the self-energy of the disorder in a coherent state path integral on the Keldysh time contour [7]-[10j. 
This equation finally results after Hubbard-Stratonovich transformation (HST) of the repulsive interaction 
and integration over the remaining bilinear condensate fields il'*{x, t) . . . ip{x, t) [HI [32] ■ We introduce a 
'Nambu'-doubling of the fields for the HST's so that the anomalous terms are strictly incorporated into the 
selfenergies and the saddle point equations [121 IE] • The self-energy of the pure interaction (as the direct 
product of M times sp{2) and as subalgebra of the disorder-self-energy sp(4M)) can be absorbed by a shift 
into the 'larger' self-energy sp(4M) of the ensemble averaged random potential ('M' denotes the number of 
discrete time intervals). This disorder-self-energy sp(4M) is represented by a larger group and also follows 
from a HST, but of a non-hermitian 'interaction' term resulting from the average over the random disorder 
[TT| [32] . After a Gaussian integral of the self-energy of the repulsive interaction, the self-energy of the 
ensemble averaged disorder only remains as a matrix valued field depending on the two branches of the 
time contour. The exact saddle point equation, obtained by a first order variation, therefore only consists 
of the self-energy of the disorder. Apart from the purely bosonic constituents of the disorder models the 
presented approach is reminiscent of the Bogoliubov-de Gennes equations and the Nambu-Gorkov Green 
function formalism in the theory of superconductivity [14j-|16). Solutions of the saddle point equation can 
be achieved by continued fraction of the disorder-self-energy starting from the free Green function [ITl [18] . 
The iteration process for the continued fractions can be simplified in the presence of spatial symmetries or 
time independence. It is the imaginary increment —i Sp {s± = ±e, e > 0) in the Green function of the 
disorder-self-energy which determines the solutions of the self-energy of the saddle point equation. In fact 
the original coherent state path integral, consisting only of the bosonic fields for the atoms, is only defined 
by introducing this imaginary increment —i Sp on the time contour tp. This imaginary increment can 
be considered as a kind of regularization and is necessary even without zero eigenvalues in the exponents 
of the time development operators. This follows because Gaussian like integrations (jl.ip as part of a 
time development operator only lead to the absolute values of eigenvalues with respect to the sign of the 
considered time intervals. Therefore, the important sign or the information about bounded and unbounded 
states would be lost under time reversal without the infinitesimal imaginary part —i Sp 



J d[tps{tp)] exp I - J dtp ^ V^K^p) (^-^£p + l^'3^ i'sih 
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Sp^±=±e, (e>0); Aip=± = ±At, (At > 0) 



The imaginary increment has therefore also to be taken into account in coherent state path integrals without 
disorder in order to distinguish between advanced and retarded Green functions which can appear after 
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transformations and integrations over field variables. One can improve the expansion of the actions in the 
final path integral to quadratic and higher order in the single self-energy matrix of the disorder for the 
inclusion of fluctuation properties around the saddle point solution. 

In section 2 we describe the average of a coherent state path integral on a time contour [19] and consider 
the two cases of static and dynamic disorder. In section 3 the Hubbard Stratonovich transformations are 
given and the exact saddle point equations are derived for static and dynamic disorder. We outline the 
various steps for solving the saddle point equation via continued fractions and indicate the reminiscence 
to many-body theory with anomalous terms [S]- [IS], [20] -[23]. Section 4 lists the relevant observables 
obtained by differentiation of the generating function with respect to a source term. Section 5 contains 
a summary and points out the extension of the mean field approach in sections |2]|4] for d = 3 spatial 
dimensions to the two dimensional case. Since the mean field approach is less applicable in the case of 
d = 2, a nonlinear sigma model can be preferred, following from spontaneous symmetry breaking and a 
gradient expansion of a determinant [24] with additional 'Nambu'-doubled pair condensate terms in the 
self-energy. However, in the d — 2 case different kinds of HST's have to be performed for a corresponding 
nonlinear sigma model as for the d = 3 mean field approach (see section |5|) . The saddle point equations 
in d = 3 can lead to nonanalytic behaviour because of the singularity of the density of states at vanishing 
wavevector. In the spatial d — 2 case one has to apply the Weyl unitary trick for the parametrization of 
the disorder-self-energy into densities and anomalous terms with a decomposition into a subgroup U{2M) 
and coset part Sp{4:M)\U{2M) so that a nonlinear sigma model can be achieved after a suitable HST and 
gradient expansion [25] -[27], [28]. In the strictly one dimensional case d = 1 one has to take into account 
large fluctuations so that the transfer matrix approach of disorder-ensemble-averaged generating functions 
can be chosen as the prevailing computational method for density and correlation functions |29j-|31j.[19|. 



2 Coherent state path integral 

2.1 Averaging methods for zero and finite temperatures 

The Hamilton operators for the disordered bosonic systems with the Bose field operators {[tps, i's'] ~ ^s,x'] 
i^x, V's ) contain the trap potential u{x), the quartic, repulsive contact interaction with parameter Vq > 
and the kinetic energy term with mass m. Two kinds of random potentials Vi{x) (static disorder) and 
Vii{x,t) (dynamic disorder) are introduced separately and result in the two different Hamilton operators 
Hi{'tpg ^i^xiVi) (|2.ip and Hiiii'i ,i}s,Vii) (|2.2p so that we consider two models I and II in parallel. We 
examine these two models of disorder at zero temperature with coherent state path integrals and incorporate 
a chemical potential or reference energy ^o- Furthermore, a U{1) symmetry breaking, hermitian source 
term with jip-sit) is included for the creation of a coherent BE-condensate wavefunction. Since we also 
consider pair condensates of bosonic atoms on the coset part Sp{4:M)\U{2M) of the disorder-self-energy, 
we have to define a hermitian source term j-ip^p-sit) for the creation of anomalous terms as {^g{t'^) tpsitp)) 
and {i^g{tp) i^^itq)) (for the contour times tp,t'^ see relations (|2. 1012. lip ) 

Hiii^iJx, Vi)^J2'^i(- + "(^) - + '^^(^) + ^0 i'ii'x]i'x+ (2.1) 
+ {^i-A*") + U;sA}) + \Y yH^A^) ^^^^ + ^ti't U^At) 



HiSt,i^s, yii)=Y.^ti- + - A^o + Vii{x., t) + Vo ^+4-") 4+ (2.2) 

X 

+ Y {^i'-d*) + UAt)) +\Y1 yh-d^) ^s^s + ^ti't H^At) 
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Especially, in the case of disordered systems, the second moments of the Gaussian probability distributions 
of the random potentials Vi{x) (I2.4|l and Vii{x,t) (12. 5p have to be normalized in such a manner that the 
broadcnings of the eigenvalue spectra of the hermitian, ordered parts in Hj{iil;i ,^pg, Vj), Hjj^ipi ,^ps, Vji) 
p.l|2.2|) remain finite under the corresponding ensemble averages in model I and II. Therefore, the spatial 
integrals J^s (|2.1l2.2p are given in a volume normalized kind as in (|2.3p for a system size of L"^, and the 
energy values or frequencies of parameters and fields are related by the scale of the inverse time interval 
r2 = 1/At. We have to scale the actions in the generating functions at appropriate steps for the derivation 
of the saddle point equations in model I and II during a considered time development between times 
—To/2 < t < +Tq/2 (|2.3p . In the cases of model I and II with static and dynamic disorder, the suitable 
normalizations of second moments for Vj (x) , Vn (x, t) of Gaussian white- noise distributions are obtained 
by the relations (|2.4p and (|2.5p with the parameters Rj and Ru, respectively (dim[i?/] — [energy • time]; 
dim[i?jj] = [(energy)^ • time]). The energy scale hn, follows from the discrete time steps, and the parameter 
Afx is the total number of discrete points in the d-dimensional coordinate space of volume L'^ with discrete 
spatial intervals Ax 

?2 02 



.. .. ^ L 



Vi{xi) Vi{x2) = fei,x2 ; Afx = ( ^ ) static disorder (2.4) 



Vii{xi,ti)Vii{x2^t2) — Rji (J^i.Sa '^(^1—^2) dynamic disorder . (2.5) 

The generating functions Z[J ^ Vi] , Z[J ^ Vu] of the disordered systems in model I and II at zero temperature 
are represented by coherent state path integrals (|2.6|2.7p of the unitary time development operators p.8|2.9p 
with Hi{ipi,ips,Vi) dH]) and Hiiiipi ,i^s,Vii) on the Keldysh time contour ((TTCTl) [32]-[SH|. We 

have to take account of the negative sign in the backward propagation of the Keldysh time contour (|2.10p 
and therefore introduce the metric {rjp — p, p = ±) (|2.1ip of the contour time tp = t±. This metric 
iVp = Pj P — ^) (|2.1ip will frequently occur in the remainder because the ensemble averages of the 
disordered systems I, II couple the two '±' branches of the contour time in contrast to ordered systems 
where the self-energy fields depend only on time arguments with a single branch of the contour, respectively 

m 

Z[J,Vi] = {0\Uii^To/2,+To/2-Vi;J) Ui{+Ta/2, -To/2;Vi; J)\0) (2.6) 
Z[J,Vii] - {0\Uii{-To/2,+To/2;Vii;J) Uii{+To/2, -To/2;Vii; J)\0) (2.7) 

ft 



Uiit,~To/2;Vi;J) = T cxp 1^- J ^ ^^dr Hiii^+ Js,Vi; J)"^ (2.8) 
Uii{t,~To/2;Vii;J) = T exp 1^^ ^ J ^ ^^dr Hii{iP+ Js,Vii; J)^ (2.9) 

dt+...+ dt = / dt+...- dt (2.10) 

-00 J +00 J —00 J —00 

/+00 

dtp r]p...; rjp = <{ 7y+ ^ +1 ; 77, = -1 j> . (2.11) 



The coherent state path integrals Z[J',Vi], Z[J',Vii] (|2.6l2.7p at zero temperature are normalized in 
the case of vanishing 'exterior' source term which allows to obtain observables from differentiating 
Z[J, Vf], Z[J, Vii] (|2.6I2.7P by J (compare sectiongl). The property of normalization ofZ[J, Vi], Z[J, Vn] 
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(|2.12|2.13|) is guaranteed by the unitary time development [//, Uu (|2.8|2.9p in forward and backward 
H- ' direction on the time contour (|2. 1012. lip . The presence of the source fields j,p-x{tp) and j^^-x{tp) creates 
Bose particles from the vacuum states |0), (0| with the corresponding coherent state fields i^sitp), i^gitp) 
and nonvanishing anomalous terms {ip3{tp) ip3{tp)), {'ipsi'^p) i'si^p))- However, we have also to require in 
final relations for observables that the source terms j-^-^sitp), j-4>ti,;sitp) have the same values on the two 
branches of the time contour. This is defined by relation (|2.14p with the vertical line and has to be added to 
the generating functions. Therefore, the required normalization property of ensemble averaged disordered 
systems is fulfilled and possible problems with limits in replica field theories or super-symmetric extensions 
are circumvented [3S1[3H]- According to the property of normalization at zero temperature (|2.12l2.13p . the 
Gaussian ensemble averages in model I and II ()2.15|2.16p are well defined and can be transferred to other 
physical problems with disordered parts for generalized coherent states (as e.g. S'J7(2)-coherent states [35] ) 



Z[J = 0, Vi] 
Z[J = 0, Vii] 



= 1 
= 1 



(2.12) 
(2.13) 
(2.14) 



Zi[J] 



(0|?77(-To/2, +ro/2; Vf, J) Ui{+To/2, -To/2; Vi; J)\0) 



(0|[///(-To/2, +ro/2; Vn; J) Un{+To/2, -To/2; Vn; J)\0) 



(2.15) 
(2.16) 



The normalized unitary time development at zero temperature (|2.6II2.16|) has to be modified in the case of 
a finite temperature. We briefly describe the suitably normalized generating function for finite temperature 
in the case of model I (static disorder, compare with the Hamilton and unitary time development operators 
()2.1l2.8p for model I). The inclusion of the grand canonical statistical operator exp{—/3{hi — /i N)} with 
hi{'ipi ,ip3,Vi) (|2.18p as part of Hi{iil:i ,ipg,Vi) ()2.17l2.ip and its appearance in the denominator with 
the trace ZfslVj] ()2.20p of the total generating function Z[J',(3, Vi] (j2.19p lead to an expansion of a large 
(n ^ oo, 71 S No > 0) limit with Zp\Vi\ ()2. 21112.23]) . This follows from the representation of the inverse 
of Zf}[Vi] (|2.2ip by an exponential integral with auxiliary integration variable x and the Taylor expansion 
(|2.2ip of the exponential exp{— x Zp\ViW the integrand with x G [0, oo) 



Hiii^tAs. Vi) = hi{ij+,ijs. Vi) + ^ {j^it) i>s + J4,;s{t)) + 

X 



(2.17) 



hi{'^t,^s,Vi) = Y.^4(^-^^ + u{x)- tio + Vi{x) + Voi^ti^s^^3 (2.18) 



Tr 



Z[J,f3,Vi] 



exp{-/3(;ij - HN)} Uii-To/2, +To/2; Vi; J) Uii+To/2, -To/2; Vi; J) 



Zf3[Vl] 



(2.19) 



Zp[Vi] = Tr[cxp{- p{hi - ^iN)} 
— — -r = / dx exp{-x Zp[Vi]} 

Zl3[Vl\ Jo 



Z[J = 0, 13, Vi] 



(2.20) 
(2.21) 
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The ensemble averaged generating function Zi[J,l3\ (j2.22|l of Z [J, /3, Vi] (I2.19|) for finite temperatures 
is given by the sum of averaged generating functions Z/^„ (|2.23p with increasing number of fields 
in so that the large (n — > oo, n € Nq > 0) limit of field theories has to be considered. The 

unitary time development operator [// is determined by the relation (I2.8|l with the operator Hj^ipi , ^pg, Vj) 
(12. 1712. ip which includes the symmetry breaking source terms and additionally the 'exterior' source variable 
J' for the observables 



Zi[J,f3] 



i-xY 



Zi,n[J,l3\ = Tr cM~P{hi- ^iN)}UI{-To/2,+TQ|2■VI■J) 



(2.22) 
(2.23) 



X Ui{+To/2,-m2;Vi-J) 



Tr 



exp{~/3(/i/ - ^iN)} 



However, we restrict in this paper to zero temperature for the models I, II. A coset decomposition into 
densities and pair condensate terms is preferable in d = 2 spatial dimensions with a gradient expansion 
of a determinant leading to the Goldstone modes of a spontaneous symmetry breaking within a nonlinear 
Sigma model [13 US]. 



2.2 Ensemble averages in model I and II 

We list in relations (|2.24|2.25p the coherent state path integral representation of the unitary time devel- 
opment operators Uj, Uu (|2.8l2.9p in ZlJ^Vj], Z[J,Vii] (|2.6l2.7p at zero temperature for the disorder 
models I and II 



Z[J.Vi] = j d[V's(ip)] 



exp 



2h 



X exp j - I / dip ^ ij't{tp) hp{tp) + Viix) + Vo V-K^p) ^sitp) ^x{tp) 



X exp 



X exp 



h 
I 

2h 



X 

dipE 

jipip^xi^p) 



(2.24) 



Z[J.Vn]^ j d[#(ip)] expj-^l AtpAt'^Y^^f{t'^) JI^A^'a^tp) n{tp) 

^ x.x' 

X exp <^ - / dip ^ ^*x{ip) Ki^p) + ^iii^^ t) + i^sih) i^sitp) tpx{tp) 



(2.25) 



X exp ■{ — I dt. 



pE [■?'V';s(*p) ^^(*p) + V's(ip) U-M^p) 

X 

^ / dipE [ih;Sp) ^dh) V'x(ip) + #(ip) #(ip) jvv;s(*p) 



X exp 



The one-particle parts are given by hp{tp) + Vi{x), hp{tp) + Vii{x, t) which consist of the common Hamilton 
operator hp{tp) (|2.26l2.27p with the kinetic energy, the trap potential u{x), the chemical potential /io, the 
energy operator —ih d/dtp of the corresponding branch of the contour time and the random potentials 
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Vi{x), Vii{x,t). Note the inclusion of the small imaginary energy increment —t £p (|2.26p on both branches 
of the time contour which allows the selection between advanced and retarded Green functions. This 
imaginary increment —i Sp determines a direction for the time development so that the coherent state 
path integrals ()2.24|2.25p are well defined for all kinds of effective energies whether vanishing, bounded or 
unbounded 



-ih 



d_ 



— A + M a; 
2 m 



Mo 



= Sp,q Vp ^(tp - t'q) ^s,x' f>'p{tp) ; ep = r]pe ; (e > 0; 7?± ±1) 



(2.26) 

(2.27) 
(2.28) 



The symmetry breaking source terms jii!-^x{t±)i j-ipip; 

s{t±) (|2.28l2.14p for the creation of a coherent BE- 
condensate wavefunction and pair condensates have to be set to the same values on the two branches of the 
contour time in the final relations for observables with vanishing 'exterior' source J^g?'g{t'q, tp). We perform 
a 'Nambu'-doubling of the coherent state fields tpx{tp) on the time contour with its complex conjugated 
fields ^^{tp) in order to obtain also anomalous terms as (ipxitp) 4's{tp)) by a single differentiation with 
Js^^sit'q^tp). We denote this 'Nambu'-doubled field by (ip(=±)) (|2.29l2.30p (with capital '^f' instead 

of the lower-case letter '■)/)') and introduce the additional indices (a, 6 1, 2) for referring to ipsitp) (a = 1) 
or to the complex conjugated part ipg{tp) {a ~ 2). There are two possible orders of the four component 
'Nambu'-doubled field ^'^'•~^^^^(ip(^-|-)). In the listing (|2.29p one gives priority with respect to the two 
branches of the contour time so that the first two components of ^'|(ip) arc on the plus branch of 
whereas the listing ()2.30p prefers ordering with respect to the anomalous pair condensates 



contour time ordering 



^a(=l/2). s 



(2.29) 



ordering for anomalous terms 



ipsitp) (a ^ 1) 
rSp) (a = 2) 



■tjjsitp) (a = 1) 
V4(ip) (a = 2) 





(«- 


1) 


\ 




(«- 


2) 


V^j(i-) 


(« = 


1) 




V m-) 


(« = 


2) 


/ 


( Mt+) 


(«- 


1) 


\ 


Mt-) 


(« = 


1) 






(«- 


2) 




V rxit-) 


(« = 


2) 


/ 



(2.30) 



In the remainder we partially use a different notation for the coherent state field variables V's(^p(=±))j 
ipgi{t'q(=±)) in the generating functions of models I, II in order to emphasize the complete independence of 
the fields ips{t+) := '05,+ (t), ipsit^) := i ps-jt) an d also ipg,it+) := V'5',+ (i), V'£'(i-) "0^',- W concerning 
the two branches of the time contour (|2.31|2.32p . However, if classical approximations are implemented 
in the coherent state path integrals of disorder models I and II, the fields iljs{tp{=±)), ^'t,{tq(=±)) may 
take exactly the same values on both contour time branches ^s{t+)='(IJx{t-), ^/'i/(i+)=V'j'(^-) (for classical 
approximations following from variations with respect to contour fields) 



contour time ordering 

\tp{=±)) 



(2.31) 



^a( = l/2) 
x,p{=±) 



( V'5-(<+) (a - 

Ipxit-) 



tpx{tp) (a = 1) 

rdtp) (a = 2) 



1) 

2) 

1) 

2) J 



( i'3,- 



it) 



V'x,p(i) (a = 1) 
rsjt) {a = 2) 

(a=l) \ 
(« = 2) 



^Ps.-it) (a = l) 

V (a = 2) y 
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ordering for anomalous terms 

,a(=l/2) 



(2.32) 



f Tpxitp) {a = 




( ipx,p{t) {a = 






2) ) 


'[njt) (a = 





/ Mt+) (a = 1) \ / i's.+it) (« = 1) \ 



Mt-) (a = l) 
rdt+) (a = 2) 



V^s,_(i) (a = l) 
^i,(t) (a -2) 



Therefore, one can also rewrite the 'Nambu'-doubled bilinear term (|2.33p with a matrix Mg? g(t'^,tp) := 
M|f^.-p(i', t) and its integrations over contour times t'^, tp as follows (by using the metric T]p, rjq (|2. 1012. lip ) 



(2.33) 



x,x' a,b=l,2 

/ + CO 

'°° x,x' a,b=l,2p,q=± 

However, the fields ^'^(i- = +To/2) and ^'j(t+ = +ro/2) have to approach exactly the same values at the 
time +To/2 for continuity reasons. This must hold in both kinds of expressions with coherent state fields 
in (|2.33p (compare (|2.29l2.3ip for contour time ordering and relations (|2.30l2.32p with prevailing order for 
the anomalous parts). 

The Gaussian ensemble averages of Z[J', Vj], Z\J , V//] (|2.24l2.25p with second moments (|2.4l2.5p result 
in relations Zi\J\ Zii\J\ (|2.34l2.35p where the 'Nambu'-doubling has also been included for the symmetry 
breaking source terms with J°;,.g{tp), J^°g{tp) and also -J^^.git'q.tp) 

Zi[J] = f d[M-tp)] exp(-| / dtpY,i4(tp) Kih) Mtp)] (2-34) 



X exp 



2h 



X exp 



X exp 



X exp 



2h 



dipE [Jt^Sp) "^iitp) + ^^pitp) j;.4tp 

X 
X 

expj-^l dtpdt'^'^^fit'^) ji^AW mtp)] 

^ x.x' 



Zii[J] = / d[V'j(tp)] exp<^ - 



X exp 



X 

d^p E [^t-si^p) mtp) + n^ih) j^tp 

X 

exp I - ^ ^ dtp J2 "fpitp) J'Axitp) mtp) 



I 

2h 



(2.35) 



8 



exp 



X exp 



X exp 



^ x,x' 

dtpj2^o {mtp))^ {Mtp)) 

X 

#7 d^X^ ^ [mp)vpMtp)) {rsm)vgMt,) 



X p,q=± 



The 'Nanibu'-doubling J!^^g^^^\tp) (|2.36|2.37p of the source term j-,p-x{tp) is obtained in a similar manner 
as that of the coherent state field ipsitp) (|2.29l2.30p . The 'Nambu'-doubled source field for j^-x{tp) is also 
defined by its capital letter J^. gjtp) with the additional index (a = 1, 2) and can also be listed in the two 
manners of ordering (|2.29|2.30jl as the doubled coherent state field ^'^(ip) 



M = l/2)f ^ 



'contour time ordering' 



ji:;xitp) (a 

jl;xitp) {a 



1 H;x{^+) 


[a 






{a 


= 2) 




{a 


= 1) 




(a 


=2) ; 



(2.36) 



'ordering for anomalous terms' 



(2.37) 



^a(=l/2), ^ 



/ j4r,x{t+) 


(« = 


1) 


\ 




(« = 


1) 




(« = 


2) 






(« = 


2) 


/ 



The 'Nambu'-doubling of the source term j^^^p-sitp) for the pair condensate terms yields a matrix J^^.j(ip) 
with local contour time dependence which can also be ordered in the two analogous kinds as the U{1) 
source term J^.^p) (|2.36|2.37p or the field *^(tp) p.29l2.30p 

'contour time ordering' : 



(2.38) 



JThjbixi^p) ^ 



I 



\ 



jrpjp-x{t+) 
jM;xit+) 



j^lp;xi^-) 



•^ipTp;x{^p) ^ 



ordering for anomalous terms' : 

/ 



jU,xit+) 
V j;^.4t_) 






jipipixij"—^ 







jil>tp:x(j^— ) 






(2.39) 



The source terms J^.-(t±) and J'^"^.^{t±) have to be set to equivalent values concerning the two branches 
of the contour time t± in the final relations for the observables or the saddle point equation (|2.14p . In the 
remainder the equivalent notations (|2. 40112. 4^ for J^.j(ip) and J^^.g{tp) will also temporarily occur as the 
equivalent notations for the coherent state fields and 'Nambu'-doubled matrices (|2.33p 

'contour time ordering' : (2.40) 

/ (a 

j-4,:x,p{t) (a = 1) 



ja(=l/2) 

ijj;x 



{=±)) 



yj;x,p{—±) ^ > 



j;;x,p(t) (a = 2) 



= 1) 

(« = 2) 
ji>;x-{t) (a = 1) 

V (a = 2) / 
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a(=l/2). 
"^ip-x V''P(= 



'ordering for anomalous terms' 



ji>;x.p{t) (a = 1) 

(a = 2) 



.U;x~{t) {a 



2) 

2) ; 



(2.41) 



'contour time ordering' 
/ 



Jljj4>;x{^p) 



Jt/Jtll-X.pi^) 



3ipTp:x.+ {^) 








jrPi);x-it) 

ju,x,-it) ; 



(2.42) 



Jlp4>;xi^p) 



'ordering for anomalous terms' : 
/ 

V j;4r,x.Jt) 



(2.43) 



(t) 

jipii!-x.-{i) 


/ 



One has to apply the matrix form J^^.g, g{tg,tp) or its corresponding notation symbol J^^.^, q.3p{t',t) 
(|2.44p in the case of two time contour integrations (|2.45p with the bilinear fields 'i>ti^{t'g) . . . ^'^(ip) as with 
the matrix M^? Sq^tp) — g (t\t) in 



•-^t/jlp;x' ,xi^qT^p) ~ •^lptp;x' ,q;x,pi^ ~ ^P,q Vp ^i^p ^q) ^x.x' J^tp-^xi^p) 

= Sp,q Vp H^p ~ tq) ^x,x' J^il,-,S,pi^) 

j dtpAt'^Y. E J'A^'AW mtp) = 

*^ x,x' aM=l,2 

/ + 00 

'°° x,x' a,b=l,2p,q=± 

= [ ^hY. "^t'itp) j'^Sp) mtp) ■ 



(2.44) 



(2.45) 



X a, 6=1, 2 



In this section we have achieved the ensemble averages of the disorder models I, II with the coherent state 
path integrals Zi[j], Zn[j] (|2.34|2.35p . We have described the various forms and equivalent notations 
concerning the coherent state fields and matrices (|2. 29112.33]) on the Keldysh time contour tp=±. The 
'Nambu'-doubling of source fields and matrices has also been incorporated for the creation of a coherent 
BE-wavefunction and pair condensate terms 



3 Hubbard-Stratonovich transformations for the repulsive and 
ensemble-averaged interactions in model I and II 

3.1 Hubbard-Stratonovich transformation for repulsive interactions in model 
I and II 

The repulsive interaction with parameter Vq > is a common part of the two disorder models I and II. 
Its Hubbard-Stratonovich transformation (HST) to a density matrix fg^{tp) p.ll3.2p is accomplished by a 
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dyadic product of the fields in the repulsive interaction term where we already insert the 'Nambu'-doubled 
form of the dyadic products with ff{tp) = ^'g(ip) (8) 



(3.1) 



dip E i^sitp)Y iMtp)f = \j dtp X] (^Kip) Mh) + Mtp) mp) 

- \l ^tpY,n{tp)®n\ip) '^'Sp)®n^{tp)^\l ^hY,^Arf{tp)rl\h) 



We introduce the self-energy matrix cri^{tp) (13. 3p for the repulsive interaction term with Vq > in the 
models I, II with static and dynamic disorder according to the symmetries of the resulting density matrix 
rfitp) in the trace 'tr^.b' over 'Nambu'-indices (HH) 



rfitp) = 
^f{ip) = 



rSp) rSp) rSp) Mh) 

4(*p) ^i'(ip) 



bsitp) e R cs{tp) e C 



(3.2) 
(3.3) 



This self-energy matrix (Jg''{tp) (j3.3p has a local dependence in the spatial coordinates and also regarding 
the contour time. Therefore, the HST of the repulsive interaction or of its density matrix form with the 
trace ^tia.b over the 'Nambu'-indices a, b — 1,2 yields the relation ()3.4|) with hermitian action in a Gaussian 
factor of the self-energy and the hermitian coupling between density matrix f'l^(tp) and ai''{tp) 



exp 



^ dtp J2 Vo {mp)f (Mtp))'^ = exp I - ^ ^0 ^ dtp tr [ffitp) ff{tp] 



(3.4) 



d[af{tp)] exp I I dtpY^j\^f{tp)^f{tp)\ \ exp{^/ dtpY^j}_^f{tp)rf{h 

The real field h3{tp) = cr^{tp) = a'^{tp) in p.3p describes the density term of the self-energy for the 
repulsive interaction whereas the complex field cs{tp) = cri^(tp) and its complex conjugate Cj(tp) = cr'i^{tp) 
in p.3p determine the anomalous terms of the interaction with Vq > 0. 

3.2 Hubbard-Stratonovich transformation for the disorder term in model I 
and derivation of the mean field equations with the disorder-self-energy 

The HST of the quartic, non-hermitian 'interaction term' for the disorder in model I is obtained by the 
dyadic product of fields in a similar manner as in section [5TT1 The 'Nambu'-doubling with ^^titp) (g) 4'i''(t^) 
leads to the disorder-density matrix Ri''{tp,t'^) with inclusion of the anomalous terms 

J dtpdt'^Y'f'titp^Mtp) dtpdt'^Y.['^t''itp)n{tp)) =(3.5) 

^ X ^ X 

= - j dtp dt'^ ® ^J'^^?) *'-(*^) ® ) = I / dip K E ^^"(^9' ^p) 



dtdt'Y Tr 

p.q;a^h 



-^x-.pqi^J ^ ) Vq -^x;qp{^ J ^) Vp 



We explicitly list the spatially local density matrix Rg''{tp, t'^) in (|3.6l3.5p with 'Nambu'-indices (a. 5=1,2) 
and nonlocal time contour dependence with tp, t'^ {p,q — ±) according to the contour ordering of fields 
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as for ^si^p) = ^spi^) P-8|2.29|2.3ip . Note that two different notations for the disorder-density matrix 
P.6I3.7|) can be used as in the cases of ^'^(tp) = *?,p(<) (|2.29|2.31|) . J%3{tp) = J%s^p{t) (|2.36|2.40p or 
Jl%;S'd^'q^'tp) = J^''i,-S',q;sjt'^i) (12.4412.4512.3812.4211 . Therefore, we have also added in relation ^ the 
last line with the trace 'Trp_g;a,6' and the metric r/p, rjq for the disorder-density matrix in notation p.7p . This 
clarifies the symmetry relations between the matrix elements of R^^{tp,t'y) = R^'^pq{t,t') (|3.6l3.7p . In the 



listings p.6l3.7p the contour time ordering prevails for the disorder-density matrices Rf'{tp, t'^) 



x;pq 



which follow by the dyadic products ^'|(ip) 
respect to the contour time order 



n it'q)^ n,p{t)®n'q^t') ofseids 



also applied with 



Rg {tp, t ) 



-^x-.pqi^J O 



V rs-) rs'+) 



m+) rs'- 
rs~) rs'- 



-■±)) = 



r^A^i^xAn nA^rsA^') 

( MA (a = 1) \ 

rA+) (« = 2) 

MA (a = 1) 
V AAA (a = 2) J 



_ ( Aitp) {a = l)\ _ 



AAA (« = 2) 



AA) A A) \ 
A A) A A) 
Mt-) A A) 
AA-)AA) J 

A At) A At') \ 
A At) A At') 
A.-it) A At') 
As At) A At') I 



'contour time 
ordering' 



(3.6) 



(3.7) 



(3.8) 



^ts;pqit,t') (13.913.1011 has to fulfiU 



The corresponding 'Nambu'-doubled disorder-self-energy T,A{tp,A 
the equivalent symmetry relations between its matrix elements as the disorder-density matrix Ri^'Aitg), 
RA„{t,t') (I3.6|3.7p . We also consider in p.9l3.10p the two equivalent notations as for the disorder-density 



matrix (|3.6l3.7p . The disorder-self-energy (|3.9l3.10p has a nonlocal dependence with respect to the contour 
times tp, t'q and contains density related terms labeled with the capital letter 'i?' and pair condensate 
terms labeled with the capital letter 'C". The basic matrices of the disorder self-energy Ej*-(tp,<^) are 
Bg{t+,t'^), Bg{t^,t'_) and BA+Tt'^) for density related terms. The basic matrices for the anomalous 
parts are given by C'j(t+, t'^), Cs{t^,t'_) and Cg{t^,t'_). Taking into account the symmetries of R'^itp, t'^), 
one has to place the basic matrices BA+it'^), BA-it'^) and BA+it'_) for densities and the basic 
matrices C'j(t+, t'^), Cs{t-, t'_), Cs{t+, t'_) for the pair condensates in the disorder-self-energy as in relation 
p.9l3.10p . Furthermore, one has to require the symmetry restrictions p. lip for these basic matrices of 
ilj'^Api t'q), A'^x-pqitj t')- The two equivalent notations concerning the nonlocal contour time dependence are 
also tabulated in relation (|3.1ip for the symmetries of the 'B' and 'C" matrices. These two notations clarify 
the symmetry relations with complex conjugation, transposition and hermitian conjugation (see following 

examples for complex conjugation (BAp,t'q))* = BtpA,t'), transposition (BAp,t'q)) = B^A^'^') ^ 

Bs;pq{t',t) = liA'p^tq) and hermitian conjugation [BApit'q)) = ^lpg(i'i^) = Bi{tp,tq) within the two 
equivalent notations in p. lip ) 



'^f^xitp^tq) — 



^l'^x;pq{t, t') 



I BsA,t'+) 




A) 


Bx{t+ 


,t'-) 


Cx{t+ 


A) \ 




cAt^ 


A'+) 


BIA 


A) 


AA 


A) 


B*xA 


A) 




BtAA) 


CIA 


A) 


BsA 


A) 


CxA 


A) 




V ctit^ 


A'-) 


BIA 


A) 


ctA 


A) 


BlA 


A) ) 




( Bs-+^ 


{tA) 


A-+^ 


{A) 


Bx-+ 


A A 


^3: + 


A A \ 




(A) 




A A 


9; + 


A A 


BA 


A A) 






(A) 


A 


{A) 


Bx- 


A A) 


Cx:- 


A A 






(A) 


4^+- 


(A) 


C+ 


-it A) 




-{tA J 



(3.9) 



(3.10) 
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■05 (^+) ^Jx{t'- 



Bs-++{t, t') 
Bs,— {t,t') 
Bs--+it,t') 

Cx;++{t, t') 
Cs;—{t,t') 
C^x;-\ — (ij t') 



1 JD^it-L 


' +/ J 


+ 


1 i?^. 


_j (-yf, / / 




.t' )' 

' -'j 


+ 

= 


\ ^x: 


/ 


I Bgitj^ 


,t' )' 


+ 


\ ^x: 


+ -(^,^')) 


( Cx(^+ 




T 

= 


( C'^; 






■''>; 


T 

= 


( C'^; 








T 


( 


- + it,t')) 



)) 
)) 
)) 

)) 
)) 

)) ■ 



(3.11) 



We use relations (I3.5P to p. lip for the HST of the ensemble averaged disorder term I with nonfocal time 
dependence and with the symmetries foUowing from the 'Nambu'-doubhng. A 'non-hermitian' Gaussian 
factor of the disorder-self-energy E°f'-(tp,t^) and its coupling to the disorder-density matrix Rg''{tp,t'g) 
result in place of the quartic, non-hermitian interaction of fields derived from the ensemble average with 
Viix) 



exp 



2nWx jc 

' 8 ft 



X Jc 



X exp 



dtp dt'^ ^ (V's(ip) ^s(tp)j (i'sit'q) i'sit'g 

X 

dip dtj^\ tr LR|^(tp,i^) R^^{t'^,tp 

X 

d[t%{tp, t'^)] exp I - ^ 1 dtp dt'^ [^%itp^ ^'U^',^ 

^ X 

^ f dtp dt'^Y. *^ \^t<*p^t',) RfiW 

VJ^x JC ^ <^'b '- 



(3.12) 



tp) 



I Rjn 



Substituting the terms of the HST's p.4|3.12p and the 'Nambu'-doubled symmetry breaking source terms 
into Zi [J] (|2.34p , we acquire the ensemble averaged coherent state path integral of model I with the self- 
energy d-g''{tp) of the repulsive interaction and the disorder-self-energy 'S'^K(tp,t'^). Moreover, we perform 

the 'Nambu'-doubling (j3. 14113. l7|) on the one-particle terms in Zi[J'] so that relation p.l3p for Zi[J'] 
only consists of 'Nambu'-doubled parts with Sjf-(<p,i^), af{tp), *|(ip) and J^,.Sp)^ Jf'P-A^ 

p) as well as 

'f^%',xiW 

d[EJ^(tp,i;)]d[<72''(ip)] expj-^l dtpdt'^Y^^jX^t^M) ^'UW 



Zi[J] 



(3.13) 



X / d[V'j(<p)] exp 



2h 



dtp dt'^Y.'^f{t'^)Nx 

x,x' 



''-x',x\''qJ '-P7 1 



X exp 



2h 



c 



2ftv7v; 



mtp] 



The 'Nambu'-doubled one-particle Hamiltonian 'H'gf g{t'^,tp) has to take the form as in (|3.14p for a chosen 
contour time ordering with the part hp(tp) p.lSp and its transpose hE{tp) p.l6p (compare with the two 
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notations of contour time ordering in (I2.29|2.31l2.33l2.38|2.42l I2.44I2.45|) ) 



/ h+{t+) 



h-{t-) 





dtp 




h u{x) 


- Mo 


hl{tp) = 


dtp 




h u{x) 


- Mo 




dtp 









h^{t^) J 



(3.14) 

(3.15) 
(3.16) 
(3.17) 



The contour energy operator Ep — ih d / dtp p.l7p is antisymmetric with respect to transposition whereas 
the other terms of hp(tp) are symmetric under transposition. It has been mentioned in the introduction 
that the self-energy d'i''{tp) of the repulsive interaction can be considered as a subalgebra (with being the 
direct product of M times sp{2)), concerning the 'larger' disorder-self-energy Ej!'-(tp, i^). This disorder- 
self-energy can itself be regarded as an element of the symplectic Lie- Algebra sp(4M), with respect to the 
number of independent parameters which is given by AM ■ {AM + l)/2 [Tfl [55], [551 [3Zj; (The parameter 
M e N > denotes the number of discrete time intervals or steps during time development between 
—To/2 < tpj < +Tq/2 for times of a single branch of the contour {p — fixed ± value, < j < M — 1, 
tpj — At ■ j). This important observation allows to shift the disorder-self-energy Yj'f!^{tp,t'^ by the self- 
energy (J'g'itp) of the repulsive interaction p.lSp . We can also transfer the source term J^^.g{tp) for the 
pair condensates as a subset of the coset part sp{AM)\u{2M) of the symplectic Lie- Algebra sp{AM) to the 
disorder-self-energy p.l9l3.20p . Therefore, S^^(^p,^'^), coupled to the bilinear fields 'i'p\t'g) . . . «'^(tp) in 

Zi[J^] (|3.13p . can absorb the self-energy of the repulsive interaction and the source matrix for the anomalous 
terms 



t%{tp, t'^) + 2-|-^ Sp^q Vp 5{tp - t'^) af(tp) 



^/^^(^Pi ^q) ^S,x' ~ 2 



Ri n 



jab (f fl \ 



'^I;xi^PT^q) Jlj 5p^q rip 6{tp t^j) J^^.^{tp) 



(3.18) 
(3.19) 
(3.20) 



After these shifts of SJ^(ip, t'^) in Zi[J] (jXT^ . we remove the 'Nambu'-doubled fields (i^) . . . *s(ip) 
by integration. According to the doubling of the fields, we obtain the square root of the determinant 
with the disorder-self-energy , the one-particle Hamiltonian 'H'^f g{t'g,tp) and with the 'exterior' source 

term >Jgf^g{t'g,tp) for generating observables by differentiation. The shifts of Sj^-(ip,t^) as in p. 18113. 20]) 

have eliminated the self-energy ^'^^{tp) of the repulsive interaction and the source matrix J^^ ^{tp) from 

the determinant. An additional Gaussian factor of cr'g'itp) and Gaussian couphng terms with E°f'-(ip, t^), 

j^'^.^itp) result instead of the appearance in the functional determinant. The Gaussian factors of the self- 



energy d'g''{tp) disappear completely from the coherent state path integral Zi[J'] (|3.2ip after integration as 
in p.22p so that the disorder-self-energy Ej''-(tp,t^) remains as the only integration variable 

Z;iJ] = J d[E?^(t„ t'g)] exp I - ^ dtp dt'g l\ [^^xitp, t'q) ^%it'q, tp)] } (3.21) 
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X exp 



X exp 



[ 1 VJT. 
1 AC 



2R]n j_ 



, a,b 
X p=± 



tr In 

a.b 



Mx 



+ - 



"iVTTxh 



^ a,b=l,2 



AC 



1 i?/^? 



2 vA/^ft 

d[af (tp)] exp 



x' .X 



tr 

°° X p=± ^ ■ 



■Hp 



hyJJTx 



Ri 



Jibjb;xi^p) 



Vp 



X exp 



exp 



X exp 



tr 



ba ( 



'/P ■4)ip;xV'P , 



" 2 hRin 



d^EE 



— , -I _ Rjn 

X p=± I -L « T7p A/-^2?iVo 



tr 



dipE 



tr 

a, 6 



^/fx(^P'^p) ^/tx(^P'^p) 
'^'l';xi^PT^p) J^M-^xi^p) 



X exp - 



1 AC 



2i?| rj 



+ 00 



d^EE 



X p=± [^ iTlp M^2hVo 



tr 



'^ibjb-xi^p) '^tbih-xi^p) 



(3.22) 



Finally, we obtain the ensemble averaged generating function Zj[J'] p.23p for the disorder model I at 
zero temperature. It only contains as single integration variables the disorder-self-energy matrix elements 
tq). The generating function ZjIJ] p.23p consists of Gaussian factors with Ej''-(ip, t^), one func- 
tional determinant and the bilinear source term with ■ ■ ■ J^-s{tp) for the coherent BE-condensate 
wavefunction. The final expression for Zi[J'] is listed in Eq. p.23p with the matrix M^.^ g{t'g,tp) p.25p 
as an important ingredient apart from the Gaussian factors with the complex parameter fip^^ p.24p 



Zj[J] = / d[S?^(tp,t;)] exp - — / dtp dt'^Y. ( 1 



Sp.„ 6{tp-t'^) (I) 



tr 



^/fxl^pi^n) ^/°?(^0'*p) 



r 1 

X exp < - 



t 2 hRin 



dtp ^ (1 - tr \tf^Sp, tp) 4%dtp) 

— a.b 



(3.23) 
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X exp 



1 



2Rin ^ 

J rt— 



x'x' a,&=l,2 



2 '^P 

-1 



(*) ( 

< 3?(m^^^) < 1 



(3.24) 



^^hx'.si^'q^tp) — Vq[T^S',x{iq^ip) + 



— ^P^Q ^P ^i^P ^q) ^x,x' 



\ 



Mx 2y^h 

( h+{t+) 

hl{t+) 

h_{t_) 

V h^it-) J 



•^S^xi^q'^p) , 1 Rl^ 



Nx 



<>x',x ^'l^sit'q, tp) )?7p . 



(3.25) 



In order to derive the saddle point equation, we have to perform the first order variation of the actions 
in Zj[J'] (|3.23p with respect to the disorder-self-energy E/^5(ip,ig). In principle one can continue the 
first order variations of the actions in Zj[j'] (|3.23p to second or even higher order variations as a kind of 
functional Taylor expansion with the disorder-self-energy. We restrict in model I only to solutions following 
from the first order variations (5Ej''-(tp, tj^) and have to scale all parameters and self-energy matrix fields 
to dimcnsionless values. This scaling to dimensionless quantities is given in relations ((X^ to (^37)) El 



dim[H|?,^(t;,tp)] 



'^x',xi^q^ ^p) 



[energy] 
[time] 



(3.26) 
(3.27) 



'^I'-xi^'q^ tp) 



Ri 



^hXi{t'q,h'tp;k) 



1 ^%M'qU^^^P-k) 



(3.28) 
(3.29) 



Moreover, we have to consider a kind of typical 'level spacing' 'e' p.30|3.3ip which follows from the funda- 
mental discreteness of spatial and time-like variables and fields in Zj [J'] (|3.23p . A dimensionless parameter 
Vq p.30p replaces the repulsive interaction strength Vq after scaling with a kind of 'mean level spacing' 'e' 



tilde ' ' over the self-energy, the one particle Hamilton operator or other parameters refers to the corresponding 
dimensionless, scaled quantity. 
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p.31|) known from random matrix theories. Additionally, we introduce the dimensionless quantity ^/ p.32p 
as the quotient of the second moment related disorder parameter Rj to the parameter Vq of the repulsive 
interaction so that the complex parameter /ip^^ (|3.24p in Z[J'] (|3.23p is determined by relations (|3.32l3.33p 

Vo 



Vn 



Ml 

Ri_ 

h 



(^o/e) 



Ri /Vo 



1 



1 '^i!;Xii^P,k) 



•^ipTplxi^p) 



'JiPi)-x'.,x,\''q.,hh,k) — 

= Sp,q rip S{tp^k,t'q l) 5xi,x'. J^^;xi{tp,k) 
ba 

'^tptp-Xi i^p-k) 



m 

~ba II '^i^-sS^'qp'^V-k) 
•^x'- ,XiV'q,h^P,k) = 



(3.30) 
(3.31) 

(3.32) 
(3.33) 

(3.34) 
(3.35) 

(3.36) 
(3.37) 



A scaling of the source terms 



has also to be included for the derivation of the saddle equation. We 



list the ensemble averaged generating function Zj [J'] 

Yfj°:^ [tq,i' ,tp^k) p.28p and corresponding scaled operators, fields and parameters 



in terms of the scaled disorder-self-energy 

in discrete 



space-time coordinates in relation (|3.38p . The generating function Zi[j] p.38p only consists of discrete 
sums with space-time points Xi, xj and tp,k, t'^ i {p,q — ± ; tp^k — k-Atp ] t'^i = ^ Ai^ ; fc, Z = 0, . . . , M — 1) 







Zi[J] = I dp?^ (tp,fc,t;,,)] exp ] - ^ ^ ^(1 - 4^)) tr [4'^,sStp,k) 4%;x.ih,k) \ (3.38) 

p=± tp,k Xi 



<^^P1 - 8 ^ E E (l - ^P-<i ^i^P-k^t'q.i) M^'^j tr rip ^tsSh.k^^i) ^9 ^tsSt'q.p^P.k) 

P''i=^tp.k,t'^^i Xi 



e^pU^EEE(i-4'^)^p ^"ikiip,k.tp,k) 4%.,xAtp,k) 



2Rj 



a.b 



^-p { - ^ E E E l^ih^. (^;./' *p^^) 

P=± tp.fe Xi 

'^-pI^ E E E ^A^^'".') ^^^'X^^^^'p^^) J^itp.k) 

p,q=± tp,k,tq,l' Xi,x'. 



Ml°x'.,xS:tqA^ip-k) = Vqi'Hs'.,Xiit'q,l,tp,k) + Js^,sS't'q,l,'tp,k) + 2 ^'l^i^'qA^ip-k) 3xi,x'. jVp (3.39) 



^In the remainder the symbol 'S{tp i^,t'^ j)' as in ZjIJ] 1)3. 3813. 35| I denotes the Kronecker-delta for the discrete times tp j. 
and t' ,. 

q,l 
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.fLe\ds{x,tp) . . . ^ ^77p . . .fields(xj,tp,fc) . . . . (3.40) 

p=± fc = 0,.-.,A/-l Si 

The sum ^ over discrete time contour variables in ()3.38p is defined without the metric rjp, rjq which 
has therefore to be included separately with the sum ^ p,q=± over the contour time branches. However, the 
metric rjp, rjq has to appear in the appropriate terms of Zi[^] (j3.38p where a contour integration is really 
performed with the negative sign for the propagation in the backward direction of a time development. 

The variation of Zi[J'] p.38p with respect to 5Yj'f!^_{tp,k,t'^ j) has to be accomplished with great care, 
implementing the particular discreteness of the contour time with its two branches and the sign [rjp — p, 
r]± = ±) for forward and backward propagation 

Sp,q S{tp^k,t'qi) ^^^^^ Y}'j'^g^{t'qi,tp,k) = 2 ^ Sp,q S(tp,k,t'qi) Jljj%sS'^p,k)+ (3.41) 



1 



_ Ri 
- Ri M-^^(t'q„tp^^-^RJ E E E E 

p',q'=± rp,_fc,,T^,^j, yii,y'., c,d=l,2 
J^;Xi^^P=+-k) — J^;xii^P=-,k) J'^4y,xS^P= + .k) — J^^-^Si^^P=-.k) ■ (3.42) 

A solution of the saddle point equation p.41|3.42p directly follows from continued fractions with the 
disorder-self-energy Sj!*-. (ip.fc, ;) as a mean field solution [17l [18]. The various steps towards a con- 
verging solution of the continued fractions simplify in the case of time independent source terms and trap 
potential. Assuming further spatial symmetries (isotropic or translational invariance), one can reduce the 
expenditure for solving p.41l3.42p with continued fractions. The continued fraction of p.41l3.42p leads to 
a solution under very general assumptions regardless of spatial symmetries or time independence of source 
terms and the trap potential. However, one has to fulfill the pole structure in the continued fractions, 
originally defined by the imaginary increment —lEp. 

The iteration m ^ m-\-l from Ej!"-. (m; tp^k,t'^ ;) to E"''-. (m-l- 1; tp,k,t'^ j) proceeds according to relation 
p.43p via continued fraction. One starts from the noninteracting, 'free' Green function with vanishing 
disorder-self-energy Y,f!g,[m = 0;tp_k,t'^ i) = and obtains a disorder-self-energy Ej**-. (m — l;tp^k,t'q i) 
with non-vanishing non-diagonal 'H — ' and ' — h' parts (see the notations for contour time ordering) 

(l - Sp,g 5{tp,k,t'q,i) /i*,')) -f 1; i^,,, ip.fc) = (3.43) 

= 2-^ Sp,g S{tp,k,t'gi) j'^^-3i{tp,k) - Rl Mj.b''^^{m;t'gi,tp,k) + 

Ri 

p',q'=± T^,y,T'^,^, yii,y'., C,d=l,2 

M'}^3,^s,{m;tp.k,t'g,i) = Vp{nt,x'itv,k,t'g^)^\RI^ts^^^^^ (3-44) 

Ji}\Xi^^P=+,k) = J^;xii^P=-,k) J^%p;xi{^P=+,k) = J'^-il,;xi{^P=~,k) ■ (3.45) 

The nondiagonal parts E^f*-. (m -I- l]tp^±^k,t'^^^ i) reappear at every iteration step due to the source field 
J^,-^(Tp'^fc') for the creation of a coherent BE-wavefunction. The anomalous parts (ipxitp) '4'x{ip)) fol- 
low from the source matrix J^^.^Xtp^k) ^ (for a ^ h\ a^h — \,2) which has a diagonal contour time 
dependence, but nondiagonal terms (a 7^ h) for the creation of pair condensates. 
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In the case of a time independent trap potential, the solution S"!"-. (tp^^, ;) simplifies to ^f^^i-pqi^k — 
t[) whose Fourier transform therefore takes the form E^**- (ui) with a dimensionless frequency u! and 

dimensionless spatial vector ^. In this case the disorder-self^^ergy and one-particle Hamiltonian reduce to 
2"^ 



. (w) (|3:49l) and Tif - (tu) (13.4713.481) in the matrix Mf-. - (w) (lOHll 



^p,g J??) Sa,b Sf fi l4x4 



V 



^7? 


h+(e) 

















h-(a / 



ah 



-I 



hp (0 = hp(0 

St. H 



BJ^^Jlo) Cj.+_{lo) Bj__iu) Cj.__{co) 
\ Ctjuj) CZJuj) B|-^_„(a;) J 



(3.46) 



(3.47) 



(3.48) 



(3.49) 



The Green function M - (w) is obtained from solving the generalized eigenvalue problem (13.501) with 
tion 

/ h(a 



the eigenfunction ^'-'"(wjv) and eigenvalue wat 



EE 

9=± f)=l,2 



V 





h (f) 









Mo 









Ri 



EE 

g=± b=l,2 



(3.50) 



c;^^(^^) B|^^++(^^) ct^_M bi^Juon) 

\ Ct^_{iON) B^.^Jlun) Ct__iu;N) B^.__{un) ) 



ah 



pq 



V'*-_^(wAr) 



R,b 



lon Vp 



R,a 



(3.51) 
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Su..,^',, = I m E E *f;(^^v') vp *|;(^A^) • (3.52) 

Note that the disorder-self-energy S^''- (w) in p.50|) also depends on the frequency lo and therefore 

has to coincide with the resulting eigenvalue lj^^ on the right hand-side of ()3.50|) . Since the generalized 
eigenvalue problem ()3.50|) becomes non-hcrmitian due to the iterations for the continued fraction, we 
have to introduce left and right eigenvectors 5" -'"(cj^,), -""{ujn) for the ortho- normalization (13.52|) . 

Assuming completeness of the eigenfunctions ^'"^'^^^'"'{ujn) p.53p . one can construct the non-equilibrium 

Green function M~i'°^- luj) (13.541) from the eigenfunctions '^'"^'^^^'"'{ujn) and eigenvalues ojat of (13.501) 

E ^f'^.t^'^^) " ^p^i ^"-^^ (^-^^^ 



M - (e>0). (3.54) 

The generalized eigenvalue problem (I3.50|) has to be considered at every iteration step to ^ m + 1 of 
the continued fraction, but reduces to the solution of the radial part of (|3.50p for a rotational symmetry. 
The non-diagonal parts of the disorder-self-energy E"''- {lo) are reminiscent of the Bogoliubov-de Gennes 
equations for superconductivity where the nondiagonal, anomalous parts are an important ingredient of BGS 
theory. The nondiagonal 'Nambu' parts (a 7^ 6, a, 6 = 1, 2) in p.50p correspond to such pair condensates as 
{jfj^^{uj) tp^pi^)) (in the bosonic case), and the nondiagonal contour time parts '-I — ' and ' — h' are related 
to quasiparticles created by defects and disorder (for a classification of various types of disorder see P71). 
Since both anomalous parts (pair condensates and anomalous disorder effects) are taken into account in 
p.50p . we have described the exact mean field theory with the 4x4 disorder-self-energy E"*"- (w). 

3.3 Hubbard- Stratonovich transformation for the disorder term in model II 
and derivation of the mean field equations with the disorder-self-energy of 
a local time dependence 

The various steps of the derivation for the saddle point equation p.41|3.42p in section can be conveyed 
to model II with dynamic disorder. Apart from the common repulsive interaction with strength Vq in 
both models, the 'non-hcrmitian' ensemble averaged quartic interaction depends only on a single time 
variable, but includes the averaging effect of V7/(x, t) in the two time contour indices {p,q — ±) (|2.35p . 
The resulting 'Nambu'-doubled density matrix Rg^p^it) p.55l3.56p . following from the dyadic products of 
'Nambu'-doubled Bose fields, therefore has only a single time variable, but two contour indices (p, g = ±) 
for forward and backward propagation apart from the 'Nambu' indices (a, 5=1,2) 



/ diE E = (3-55) 

= 1/ '^^E E {^*3{tp) Vp i'sitp) + i^xitp) r]p ^l{tp 

X p..q=± ^ 

X (^*s{tq) Vq Mtq) + Mtq) Vq ^*x{tq)^ 

= 7 d^E E i^i 

X p,q=± ^ ^ ^ 



\tq) Vq *l(i,; 
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1 POO 



X p,q=± 



B.t'' (t) 



(t) 



R-xipni^) ^ 



V'x,-W V-i+W V'x,-WV'x,+W V'x,-W^s_W V'x,-WV'x,-W 
VV'i.WV'i'+W V'i-W^s.+W ^i-W^CW 



(3.56) 



Consequently, the HST for the disorder term in model II can be performed as in relation (|3.57p where 
the doubled disorder-self-energy '^j'j.g.pq{t) G sp{4:) (I3.58P has to fulfill the equivalent symmetry relations 

p.59|3.60p as the density matrix ^|''„„(t) (|336l) 



exp 



2K^ 



exp 



/oo 
-oo ^ ^ „ L 

Vp R-xipqi^) Vq R-x-.qpi^) 



p,q=± 



(3.57) 



8 ft2 



d^E E \ 

a,b 



d[S??;5;p,(0] exp 



X P,9=± 

1 



8 fi// J-oo 



d^EE 

X p,9=± 



X exp 



4fi 



d^EE 



tr 



x p,q=± 



'Ip ^II;x;pq 



it) Vq R'x-qpit) 



^II;x;pqit) 



/ Bs,++it) Cs,++{t) Bs,+_{t) C3..+-{t) \ 

Ct++{t) Bs,++{i) Ct+Jt) Bt^_{t) 

B*s. + -it) Cs;+-{t) Bs:—{t) Cs:—it) 

\ Ci+_it) Bs.,+-{t) Ct_Jt) B,,^^{t) ) 



(3.58) 

(3.59) 
(3.60) 

Inserting the two HST's (|3.57p and (|3.4p into Zii\J] (|2.35p . we obtain the ensemble averaged coherent 
state path integral Zn [J] p.6ip for dynamic disorder and also extend the 'Nambu'-doubling to the source 
terms and one- particle part Ti^? j(i^, tp) (13.14113.1711 



(B3,++{t) , Bs,__{t) e r) (B3,+-{t) , B3.,-+{t) e C); (B3,-+{t) = Bt^_{t) 
{Cs,++{t) , Cs-—{t) , Cx-;+-(t) e C) [Cs,-+{t) = C7x;+-(t)) • 



Zn[J] = / d[I]??^,.p,(t)] diafitp)] exp 



I 1 



X exp 



-— y 



At 



dipE Vfitp)^fit 

_ a.h L 

X 

E k ^?t>9(^) ^1 ^H;x;qpit) 
a.h L 



(3.61) 



X / d[V'5(ip)] exp 



2h 



dtp dt'^J2'^+!>it'^)M, 

x,x' 

1 



njba /.I , N I ^x^At'q^tp) 
l^x' ,x\''q^ ^p) ' . r ' 



Afx 



5{t - t') 5s,s\5p.q Vp ( J7v.;x(ip) - <^sitp)) + 7^ Wl:x:,qp{t) 



mtp) 



21 



In analogy to p. 13113. 25)) in model I, the disorder-self-energy ^jj.s-pq{i) with local time dependence can be 
shifted by the self-energy a^''{tp) of the repulsive interaction and by the source term J'^^.g{tp) for the pair 
condensates 



^U-x-pqi^) ^ ^H;x;pqi^) ^ ^ Sp^q TJp J^^-gitp) , 



(3.62) 

(3.63) 



so that the matrix M^^^, g{t'^,tp) ((HTBS)) coupled to the bilinear Bose fields m°{t'^) . . .^^{tp) in Zu[J] 
P-64p only contains the disorder-self-energy '^jj.g.p„ 



Zu[J] = / d[E?J^,.^,(0] expj - -I^ Yl /°diE k ^"lUpqit) V, ^lUqpit) 



(3.64) 



exp j - i /" ^?7p V hflM^ tr In 

I ^ Jc a,b 

expj^^/ dtpdt'^Y.^^ J+^,,(t; 

x.x' a,fc=l,2 

'7g ( '^x',x{^q^ ^P 



) — tr J^^.g{tp) Jm-x{^p) 



Vq { '^x',xi^q7 ^p) 



Mx 



x^x 

•^S^xi^q^^p) . 1 



'"^ + ^S{t - t') Ss,s' ^n-s-qpii) ]Vp 

-\-OQ 



Mx 



(tq^tp) J-ip-xi^p) 



X / d[af{tp)] exp - - E / dt E - ^r;, — tr \af{tp) af{tp) 



Rji 2"'hVoJa,b 
\ - 7^ E ^ tr [^^(tp) (r;p ±fj,,.^pp{t) - 2 ^^^^-(ip; 

L ^11 p=±J -oo J '^.f' ^ 



M'i%s',x{W^ (3-65) 

Tvba ^^ \ , '^S^xi'^'q^'^p) 1 c^, ,/n c f^fca ^,\\ '^■S\xi'^'q > 
= Vq[nx',x[tq,tp) H ^ 2 ^ ^ ^Il-S'.qpW ^ ^1 Jf + 



+ S{t - t') 5s^3- 



5p,q IJp 



( h+{t+) 

hl{t+) 

h^it-) 

V h^{t,) ) 



+ :^Vq ^%x;qp{t) Vp 



The remaining Gaussian factors with the quadratic self-energy (y'g'itp) of the repulsive interaction and its 



(3.66) 



couphng to ^u-s-pqit), j^^-sitp) can be integrated out completely as in section [3^2 



e^^P i - At- E E k-''(ip) (,7, t'jls-.ppit) - 2 4;.j(ip) 

I ^^IIp=±J-°° ^; ^ 
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exp 



X exp 



'II p= 
1 1 



+ OC' 



tr 



/+00 2^ 



+ 00 



1 



2ri Vo 



tr 



tr 



X exp 2 / / / / 



Substituting p.66p into p.64p . we finally achieve the ensemble averaged generating function Zij[J'] p.67p 
for zero temperature with the disordcr-self-energy ^"^.-.^^(t) p.58p as the only field variable and the matrix 



Zu[J] = d[Ef.S;pM 



X exp 



1 1 



+ OC 



y: I dt5:(i-5,,4-))tr 

p.,q=± 



Vp ^II:x:pqit) Vq :x:qp{t) 



X exp<;i^^ 



2R 



II p=± 
1 1 



dt 



+ OC' 



X exp < - --^ ^ / d^E 1 - Jtpdh) j'Axitp) 

^ ^11 p=± "'-00 X ^ a,bL 

exp I - i 1 ^Vp E ^J, [MiIx',S',^ *p) 



(3.67) 



h^p 



1 i + 



< 5R(4^^') < 1 



(3.68) 



The scaling of the disorder-self-energy and the other energy parameters to dimensionless quantities is listed 
in relations 



'^x',xi^qj ^p) 



R 



II 



H-p 



'^x'-.Xi i^q,h ^P,k) 



S{t-t') llj^j.^^ Jt) 5{t'i,tk) Y.j^j.^^ Jtk) 



p2 _ I^II 



y ^0 

° hVlAfx 
^11 ^ 

1 

1 - t r/p 6/ 



,5(ife - 1[) 



(3.69) 

(3.70) 
(3.71) 
(3.72) 

(3.73) 
(3.74) 
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'^ii;xi i^p,k) 



•^^xp:x'.,Xi i^q,l'^P,k) 



hn^Mx 

Sp,q rjp ^itp.kjt'q l) Sg. g'^ J^^;Xi{^P-k) 
'^^°tp;xi i^P,k) 



•^iii^-^Xi i^P,k) 



3^x'. ,Xi i^q,h ^P,k) 



KlNx 
•^Sj.Xi (^qj^ ^P.k) 



(3.75) 
(3.76) 

(3.77) 
(3.78) 



The coherent state path integral Zii[J] (|3.67p is transformed to relation (|3.79p with discrete spatial and 
time-like variables of the dimensionless parameters and fields defined in 



Zu[J]^ d[^ri;x..^pqitk)] exp - 



1 1 



2 R2 

^11 p=± tfc Xi 



J^'4);Xi,p{^k) Jij;ip;Xi,p{^k) 



^fba 



X exp 



\W ^ EE(l-'^f.?<'0 \\vp^fl-,x,,pq{tk)Vqnix,,qp{tk) 
^11 p.q=± tk Xi 



1 1 



o E E E (i - ^lU-.Ppitk) j'As.A^^) 



X exp 



2 R2 

^11 p=± tk X, 

2 ?l V a.b 

p=± tk Xi 



a,b 

^II;x'.,Xi {tq,h'tp.k) 



(3.79) 



^-P{^ E E E ^'^ ..(^;.^ VP Ji,,^^p 

^ T) ri—-\- *. +' ^. ^' 



~jha 



M^y^^s^{tii,tp,k) = »/,(^wi?.-.(i;,/,tp,fe) + + ^ 4-,.- wux,.qp{tk)yp . (3.80) 

A functional Taylor expansion can be performed on the actions in p.79l3.80p with respect to '^^"^.-..^^(ife) 
as in section 13.21 We restrict to the vanishing of the first order variation in p.79p and so derive a mean 
field equation (|3.81l3.82p for dynamic disorder with a dependence on the disorder-self-energy as the only 
remaining field variable 



^ (l - Sp,, 4^^)) ^'iU.,,p{tk) = Sp,, 77p (l - 4''') 4%;^.,pitk) 
Rjj Kjj 

- M-^i:,M^,k,tp,^-^ E E E E 

p',q' = ± T^i.T^, Vi, ,y'., c,d=l,2 



(3.81) 



V,fc') V Ji>;yiX'^p',k') 
J4>-Xii'tp=+.k) = Jil!;xi(tp^-^k) Ji!'jp;xiiip= + ,k) = J^\,-^Xiitp=-,k) ■ (3.82) 

We also have to consider the originally introduced imaginary increment —i £p (|2.25p for time reversal 
symmetry breaking so that these convergence properties are consequently transferred to the disorder-self- 
energy in the continued fractions. The iteration m — > m -I- 1 of 5]j|.-..p^(m; t^^) in the continued fractions 
follows in analogy to section 13.21 and equations (|3.43p to (|3.45|) starting from the free Green function with 
^n;x,,pq{m = Q;U)^Q 



^ (l - 5p^, 4")) Sj",;5,;,p(m -t- 1;4-) - ^ 5p,q Vp (l - ^^i"^) 4%s^,pitk)+ 

Kjj Kjj 



(3.83) 
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- Mjl^\^{m;tg^k,tp,k)-i ^ 

p',q' = ± r^,,r[, y^i ,y'., c,d=l,2 

= fex' itp,k,t'i) + (ip.fc, il;) + ^ <5(tfe, <55.,5' _(m; i^)) r,. 



(3.84) 



(3.85) 



The nondiagonal parts ^^;p(?Ti + 1; tk) of the contour time also occur at every iteration step due 

to the source field J^.g,, {Tp\k') as in the disorder model I of section [X^ and the source matrix J'^'^g. p{tk) 
creates the anomalous terms. The solution '^jj-g.-pq{tk) of ()3. 83113. 85|) reduces to a time independent 



function E?r.- in the case of a time independent trap potential u(^) where ^ is the dimensionless spatial 



I I-Xi -pq 

vector. In order to obtain the corresponding Green function MT}'^^^, {rn\tp^kit'„ i) (|3.84l3.83p . one has 

___ ^^^^^j ' q^ 

to solve the eigenvalue problem (|3. 90113. 9T|) as in section 13.21 which results from the following relations 



for the matrix M ^'i''' ^ (uj) after Fourier transformation to dimensionless frequency uj 



ii;i,p;i',g^ ' 



n-l,pq 



Izl 



/ h+(a 




V 





h+(0 









h-(0 








h-(0 / 



ab 



(3.86) 



(3.87) 



(3.88) 



pq 



hp(0 = £p - % • m(0 - Mo 



hp(C) = hp(0 ■ 



(3.89) 



In comparison to the disorder model L one has also to compute the right and left eigenfunctions ^'^^^'"(ajAr) 

and eigenvalues ujn, but without a dependence on the eigenvalue of the disorder-self-energy - for 
stationary states 



EE 

q=± b=l,2 



I h(0 



V 





h (a 









-h(a 











ab 


1 ?A|-_^(wAr) \ 


R,b 










^f_(^w) 






V'*-_(wiv) ) 


q 



(3.90) 



P9- 



,T, H , fa / X 

(ujn) 

€,9 
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q=± 6=1,2 



ujn rjp 



1 























\ V'*-_('^JV) J 



R.a 



a6 



+ - 



R,b 



h(e) = -a^^-9^^+u(|)-Mo h^(C) = hp(0 (3.91) 

-^^^.^V = /dKlEE *J;(c^^v') • (3.92) 

^ U ^ — in ' ' 



p=± a=l,2 



We transfer the generalized eigenvalue problem p. 50113. 541) in section [3?2l to the case with dynamic disorder 



-l;afc 

the orthonormalized eigenfunctions (I3.92p and eigenvalues 



so that the Green function M - (w) (13.941) is also determined by relations (13.901) to (13.931) in terms of 

li',^,P',^' ,q 



vp *|;(^A') *|;',('^A') = '^f^? ^-'^ (3-93) 

M-'f. (to) - (r>0). (3.94) 

We also assume the completeness of the orthonormalized eigenfunctions (I3.93|3.92|) as in section 13.21 for 
model I. The eigenvalue problem (|3. 90113. 9^ is reminiscent of the Bogoliubov-de Gennes equations and the 
Nambu-Gorkov Green functions in the theory for superconductivity IIl]-[Il],[in]-[35. It has to be solved 
at every iteration step of the continued fractions, but is easier to solve as in the case of the static disorder 
because the disorder-self-energy is time independent for a trap potential having only a spatial dependence. 
The iteration procedure can be further simplified in the case of spatial symmetries. 



4 Determination of the observables from the derivative with the 
source term *q) 

The C/(l) invariant density 'lim5t_^^o+ (*i'^"^(i+) *|=^(i+ + (5i+))' of non-condensed atoms follows from 
differentiating Zi\J\ p.23p . Zii\J\ (|3.67p with respect to J'gA'"'~^ {t+,t+ + (5t+) and the appropriate 

normalization. Since a field operator 'ijjxi^p) and its hermitian conjugate ^g{tp) must not act at the same 
space time point due to the infinite delta function of the corresponding commutator at coincidence of time, 
a limit process (lim^tp^Op ■ ■ - t P — ^) has to be performed for the density terms at the same branch of the 
contour time. The corresponding relations for static and dynamic disorder are tabulated in Eqs. (|4.ip to 
^ with the matrices df,;,{tj>,t'^) = M^i%,3'[ipA) and O^-^-, (tp, t^) = Mf^,3,AipA) The 

one particle part 'Hg^g,{tp, t'^) is defined in Eqs. (j3. 14113. iTP ('c, d = 1, 2' arc 'Nambu'-indices as 'a, 6=1, 2') 

static disorder : see Zi[J] ((X^ with M'ifgg,{tp,Q ([X^ (4.1) 
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dynamic disorder 
na''^ (f f' ) 



see Zjj [J] with Mfi.3:g,{tp,t'^) Km 



(4.2) 
(4.3) 



-zhn'jv^ d^^''-'^''='(t^ + 6u,u) + j dt; dt;' ^ J] 



Xi .X2 c,c/— 1,2 



(4.4) 



^lfx,x'{^PT^'q) — Vp (^S% (^p; ^q) + 2 y/j\f~ Ji ^?fs(*P'^«' 
M}f,x,x'{tp,t'g) = ^p(H|j?'(ip,i;) + ^5(i-0 SJ/;x;p,W)'79 



^79 



J^;xi^+) — J'ip;x{^-) 



J4ith:xi^+) ^ ^M-xi^-) 



(4.5) 

(4.6) 
(4.7) 



We can disentangle the general relations (|4.1j|4.7p by Fourier transformation to energy momentum space 
for a time independent trap potential. The non-equilibrium Green functions 0^]^f'^{tp,t'^) become more 
transparent in energy momentum space O'^'^ (uj) (|4.8II4.15|) and are analogous to the 'Nambu-Gorkov' 

Green function formalism for superconductivity [2] -[S], [20] -[IS]- In the case of d = 3 spatial dimensions, 
we have to consider the zero momentum state in the summation over wave- vectors X]fc • ■ ■ explicitly (j4.8p 
[21 [3]. The mean field equations p. 41113. 45)) . p. 81113. 85)) are mainly applied in three spatial dimensions 
whereas the two dimensional case should be treated preferably by spontaneous symmetry breaking for the 
derivation of a nonlinear sigma model and requires different HST transformations than the ones described 
in this paper (see section [5] and [37|.[28|) 
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. (fields, fc 7^ 0) 
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lhn^(vl/+'^=^(i+) n=\t^+^t^))''^=^^^^xY. I ^ eM-^St+Lu} 6z^j;:^l=^^ic.) 



+ n^Tgj^^yy y / ^ ^ (c^") Or'^'b'r'K) exp{ 

X exp{z((^- fc') -x-icj'- u") U)} dz^t:(^') j;,uS<) 



A;, + ;A:i ,p 



(4.9) 
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(4.10) 
(4.11) 

(4.12) 
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hp(fc-fc') 



(4.13) 
(4.14) 



*p,q Vp ^c,d 
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h_(fc-fc') 
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h^(fc-fc') / 



cd 



pq 



2m 



Ep-Mo 5j:j:, + u{k ~ k') 



(4.15) 



If we further assume spatial independence of the source field - (cUp) and of the trap potential u = 

u{k — k' = 0) and also a constant creation rate with J"; - {uj' ) (|4.16p . the U{1) invariant density terms in 

()4.4|4.9p reduce to integrals over energy momentum space and a contribution from k = 0, uj = which is 
related to the 'negative' density of the coherent BE-wavefunction 



27r i 
tp;k^ To 



^fc,o 2 



J e R 



(4.16) 



static disorder 



do; 



(4.17) 
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Ak=0;p'q' 
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+p 



hm (*+'''=i(t+) vl/-i(t+ + ^t+)) 



dynamic disorder 



X 'I^V,^' ^P' l4x4 
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2m 
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exp{-z Ji+ X (4.18) 
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The coherent BE-wavefunction ipBEcix,t+) is obtained by differentiating Zjjj[J'] ()3.23l3.67p with respect 



to the U{1) symmetry breaking source field J^'^ ^(^+)- The general case for a spatial and time depen- 
dence of the coherent BE-wavefunction is listed in the following Eq. (|4.19p where the Green function 
Os^'^'-'it+.fp) refers to M^j^^.Ju;) |4.l0|4.12D and to M^lr.^^.j:, Joj) (|4.1l|4.l3D 



(4.19) 



We can transform the above relation (|4.19p to energy momentum space which simplifies for a homogenous 



translation invariant system with u = u(Ak = 0) and jii^p) — if ^i^) % o 2 ' ^ l|4.16p 
i;BEc{£,U)^AfMTo f E exp{z(fc.f-^'t+)} 



(4.20) 



static disorder 



Ri 



EE] V,?' V ^c',d' l4x4 ( « £p' + Mo - W ) + 

— l;a— l.c 



EEC 



(4.21) 



^ V SlLo;,','(- ^ 0) ^ 0,0; ^ 0) l^ToAA^I 

dynamic disorder I ^ r ^ / \ / a nn\ 

= 2^ 2^ S V,9' V I4X4U V + Mo - w )+ (4.22) 

-| N — l;a— l.c 

- 2 ^l£=5;P'«' (fc = 0, ^ 0) ryp f}ToAA,| 

In the thermodynamic limit j ^ 0, a finite coherent BE-wavefunction ipsEC for a homogenous system 
remains if an 'effective zero eigenvalue' appears in the denominators of (j4.21|4.22p . The order of magnitude 
of the BE-wavefunctions V'bbc (|4.21l4.22p can then be estimated as follows 



IpBEC 



n To j 



I e+ + ('effective zero eigenvalue') 



finite value 



(4.23) 



Using the properties of non-equilibrium Green functions [33 , the C/(l) invariant density of non-condensed 
atoms 'lim5f^^o+ ' ~^(*-i-) ^^^^(^-i- + '^^+))' contains the density \'ipBEc\'^ of the coherent BEG wave- 
functions (|4. 20114. 2^ which is subtracted from the total density given by the first terms in (|4.17l4.18p 



hm {^p''=\t+) ^Sr\U+St+)) 



static disorder j-r^ k r \ ^ f r ci ~i 

= iMlNx2_^ / TUT exp{-idt+cj} 



Sp',q' Vp' ^c',d' l4x4 



hU! + t Er,' ~ I — Mo + U 
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(4.24) 
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Sp',q' Vp' ^c',d' l4x4 



\k\ 

2m 



fJ'O 



2 '^P' ^Ak=d;p'q' ^1' 



_l:a=l.b=l 
(fc,C.) 

+ + 



V'. 



BEC 



The corresponding relations for the bosonic anomafous or pair condensate terms {ip3{t+) 4'x{i+ + ~ 
+ 5t+)) can be taken from Eqs. (gj]) to ((4l^ by setting the index '6 = 1' in 
these equations to the value 'b = 2'. This follows from the 'Nambu' doubling of fields where ^'i'''~^(t+) 
is equivalent to The limit process lim5t^_^o+ need not be considered for the anomalous parts 

{tpg{t+) ipx{t+ + <5<+)) whereas the limit process is of central importance for the density term ()4.4I4.9|) 
lim5t^_^o+ '4'xit+ + ^t+)) because one must not compute with the operator ips{t) and its hcrmitian 
conjugate ipgit) at the same space time point in coherent state path integrals of many body theory. 



5 Summary and conclusion for d = 2 spatial dimensions 

It has already been mentioned that the derived mean field equations ()3. 41113. 45|) . (j3. 81113. 85|) for static 
and dynamic disorder in sections 13.21 13.31 are particularly applicable for d = 3 spatial dimensions because 
of the singularity of the density of states at A: = 0. Since the mean field approach is less applicable in 
d = 2 dimensions, we briefly describe and point out an alternative method [371 US HI] which extracts the 
Goldstone modes from a spontaneous symmetry breaking in a nonlinear sigma model for the anomalous 
pair condensates. According to this procedure of spontaneous symmetry breaking, the HST transformation 
has to be considered in a different way with so-called 'hinge'-functions <^S^!^^(i) (|5.1|5.6p . 5(Tt°'{tp) (|5.13^ 
I5.17P as subgroups of the 'larger' symmetry groups of the total disorder-self-energy and the self-energy 
of the repulsive interaction. Using these modified HST's, a nonlinear sigma model can be derived from 
spontaneous symmetry breaking and a gradient expansion for the Goldstone modes. After introducing the 
block diagonal densities as 'hinge'-functions, they can be eventually removed from the determinant and the 
part for the coherent-BE wavefunction with the bilinear source field J^.g{tp)- They remain in Gaussian 
integrals and parts of the factorized invariant measure and can be eliminated by integration with a coupling 
to the source matrix for the anomalous terms. In order to acquire the HST for the disorder-self-energy of 
model II, we introduce the diagonal self-energy cr^^^^ {x, t) (jS.ip and the 'hinge' parts 5Yi]j^pq{x, t), 5Yi'^^^{x, t) 
(|5.3l5.4p and also the terms 5Y}g.^^[x,t), 5Y?g.j^^{x,t) (|5.5p with 5cg.pq[t) (|5.2p for the pair condensates 

^ ^ 'hinge' functions : (5Sypq(i), (5.1) 

5cs,pq{t) e C; 5B3,++{t), 5B3,— {t) e 5Bs,+.{t) e C (5.2) 

.yii _ f SB3.,++{t) 6B3.,+-it)\ 
o^x;pq{V - 6Bt^_{t) 6Bs,__{t) J 

"^x;pq 



"^xipq 



One can generalize from the disorder model II to model I with a double time dependence in the disorder- 
self-energy. In the case of a stationary state, the self-energy of model I with static disorder obtains a single 
energy dependence after Fourier transformation because of the reduced dependence to the difference of the 
two times. Therefore, we can replace the single time dependence of the disorder-self-energy in model II with 
dynamic disorder with the single frequency dependence of model I with static disorder under restriction 
to stationary states. The corresponding HST with 'hinge' functions is listed for the ensemble averaged, 
non-hermitian interaction of dynamic disorder in the following relation 
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dtJ2 Tr 

^ p.q-.a^h L 



X exp ■ 



4J1 



dt^ Tr 

, _ p.q-.a,o 



?7 




Note the minus sign before 'Sf^pq{t) (|5.6|5.7|5.12p and the tilde ' ' of 6J:'t'^J'^{t) (IS7|) which refers to 'anti- 
hermitian' anomafous parts '^^^^^^(t) = i ^'^'g.pqit) (a 7^ b) in this section [51 We have also included a 
second metric K (|5.6l5.9p apart from the metric K^^ = fjp 5p^q 6a.b (|5.8p for indefinite orthogonal symmetry 
concerning the two branches of the contour time. The diagonal matrix K (|5.9|) is the appropriate metric 
for the symplectic Lie algebra which changes the disorder-self-energy Sl^g^pqit) K to an element of sp(4), 
thereby fulfilling the exact commutation relations (with antihermitian coset parts !) 



m-pq{t) = 



pq 



Sa,b ^p,q Vp ^ 



S^i'^pqit) = I 5Y.%q{t) (a ^ b) 



a,b |, 



-1 ; -1 



(5.7) 
(5.8) 
(5.9) 



^s]pq{i) 



(5.10) 
(5.11) 
(5.12) 



'^rI t) Vp 5p,q + ^^Wpqit) 

-afjx,t)inp5p,q + 5t^i^p^{t) 

In a similar manner the repulsive interaction term can be transformed with a diagonal self-energy (jy^ (x, tp), 
'hinge' functions SaV-{tp) — Sa'i'^{tp) and anti-hcrmitian anomalous terms 5alg'^''\tp) = i 5a^g'^''\tp) 

. 2 / \ 2 



vo {rSp))' (Mtp) 

J d[CT(^°^(^,tp)] exp|^-^^dtp^o-[^^^(f,ip) cr|^^^(f,tp)| 
J d[5atitp) H] exp I dtp ^ [s^s\tp) « 5^1" (^p) 



(5.13) 



X exp ■ 

4'{tp) 

^¥{tp) 

SaTitp) 



2h 



dt 



pE 



tr 

a, 6 



-R|;pp(*) ^l?pp(0 



'^IHtp) SalHtp) 

s4Htp) -Tfitp) 



{x, tp) + Sal^ (tp) ; Sal' (tp) , Saf (tp) e R 



e R 



5aT{tp) 



Safitp) = I Safitp) (a ^ b) 



(5.14) 
(5.15) 
(5.16) 
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C 



(5.17) 



In analogy to relations p.l8ll3.2U)l . (I3.62l3.63p . we continue by shifts of the total disorder-self-energy with 
Sa~''{tp) (|5.18p and of the self-energy ay^{x,tp) of the repulsive interaction with (t^^^^(x, i) (|5.19p and also 
include the shift with the source matrix J^^,.g{tp) (|5.20p for the creation of the bosonic pair condensates 



AO) I 



■.■pq\ 



SJ^^tJt) - 2 Sp,, T^p 5af{tp) 



'^^ifpnlO ~ 2 5p^q rjp J!m)-x(J'p) 



(0) 



(5.18) 
(5.19) 
(5.20) 



After several transformations we finally achieve a coherent state path integral (|5.2ip which only depends on 
the anomalous terms, determined by the matrices Tp^{x, t) (|5. 25115. 27| of the coset part S'p(4)\[/(2), and on 

the self-energy (t[^^ {x, tp) of the repulsive interaction. We list as starting point for a gradient expansion the 
relation ()5.2ip where the block diagonal 'hinge' functions <5S|°p^(i) are still present in Gaussian factors which 
are to be removed by integration after a change to the corresponding invariant measure for 5'p(4)\[/(2) ® 

U{2). The diagonal self-energy cT^^^^(a;, t) has already been absorbed by a shift into ay^{x,tp) in the 
determinant and in the bilinear term with J^.g{tp) so that its remaining in a Gaussian factor has easily 

been eliminated by integration in Zji[J'] (|5.2ip 

1 



=exp<^ - 
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„ lJ-oo - ^ ^ a,6 L 



(5.21) 
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E E E E E -^^ 

' p',g'=± a',6' = l,2 p,g=± a,h=l,2 5,5" 

4Ut',')iKf^,l{x',t') d-]t{t'^,tp) f;^]-'^^'{x,t)i4^Sv) 



d'ls{t',,tp) = (5.22) 
= Sa,b 6p,g ds^s> 5{tp - t'q) (H;{tp) + (f , tp)) + (f -1 i^IKty^^ Jt'q.tp) + 
+ <5j,5, 5(i - t') [f-;-'"''{x', t') H;,{t) f;;;(x, t) - Sa.b Sp,, H;{t)) 

i5-H(f-i,f ) 
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Hp-\tp) = hpitp) ^ -ih-^ ~ I Ep - + u{x) - (5.23) 

H;-\tp) = h^{tp)^+ih—-iep~—A + u{x)~fio . (5.24) 

The matrix f^q{x, t) (I051) in O^^i ^{t'^.tp) (|5.22l5.2ip contains the pair condensates with matrices Y^q{x, t) 
dOi), Xpq{x,t) as the coset part Sp{A)\U{2) of 5*^(4) 

f;l{x,t) = (exp{-f/^'!'(x,i)})"' (5.25) 



11 



(?7pX+(x,t)77,)'' (o) 



= I . x2i ' ,22' I (5.26) 



P9 / pq 



One can extract the Goldstone modes of a spontaneous symmetry breaking Sp{A)\U{2)®U{2) in a gradient 
expansion of the operator 5'H{f-'^,f) with the matrix f^^{x,t) (IF?^ in 6%^ g{t'^,tp) The relevant 

parameter for classifying the various terms of the gradient expansion is the number of discrete space 
points. Furthermore, one obtains special properties of the coefficients multiplying the traces of the gradients 
with the matrix Tp^(x, t) va d = 2 spatial dimensions. Apart from the conformal invariance of the nonlinear 
sigma model in d = 2 POI, 1^ , the coefhcients reduce to one point functions in the spatial isotropic case 
which allow computations by saddle point approximations for the coefficients containing the self-energy 

The one dimensional case with white noise disorder can be preferably treated by transfer matrices of 
ensemble averaged generating functions because large fluctuations about mean field solutions may occur ^29, - 
|31j . One can also try to extend the transfer matrix approach to d = 2 spatial dimensions by approximating 
and restricting to the lowest momentum modes perpendicular to the transfer direction j31j . 
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